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Dans ce chapitre, on désigne par Æ le corps des réels R ou celui des com- 
plexes C, et par Æ un espace vectoriel normé de dimension finie sur X. On note 
(L(E) ,+,.,0) l'algèbre de Banach des endomorphismes de E. Ce chapitre est 
une application directe du cours sur la réduction des endomorphismes [1]. 


1 Exponentielle d’un opérateur 


5 ne < n < . a 
Le rayon de la série entière Ÿ°,,-0 %r est infini de sorte que cette série converge 
nm 
sur tout R (vers e*). Si u € L(E), la série D}, Ie convergera, et par 
: n 
conséquent la série ÿ},-9 # sera normalement convergente sur toute boule 
fermée B de centre O de £L(E). On définit ainsi un opérateur noté e" et 


appelé l’exponentielle de « : 


n>0 
Théorème 1 On a : 
(1j = fd. 
(2) Vu,ve L(E) uv = vu e"tv = etev. 
En particulier e“ est inversible et (e“)7} = e7% 


(3) Vu,ve£(E) udu=vu#ue = eu: 
(4) VueL(E) le“| < ellul, 
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Preuve : Seul (2) présente une difficulté. On a 


N N 
D . u" vP : u"vP 
ee” = lim ) — ) — | — Jim An avec An = ) ——, 
N—+00 n! p! N—+00 nlp! 
nm—= 


0 p=0 O<n<N 
OSpEN 
N k _ 
u+v uv? 
eutv — ) ur = lim Bn avec Bn = 5: Ù a 
k=0 | k=0 n+p=k M à 


Si l’on représente dans le plan les points de coordonnées entières (n,p) sur 
lesquels on somme, on constate que AN correspond à un carré et que Bx 
correspond à un triangle contenu dans ce carré. Par suite 


P +oo k 
Uu v u|| + [v 
| An — Bnl < ÿ d | | 5 (Ilu] _ I) 


k=N +1 n+p=k k=N+1 


et cette dernière somme tend vers 0 quand N tend vers +co, puisque la série 
D_n>0o 27 converge pour tout x ER. sm 


Théorème 2 Soit ue L(E). L'application 
f: (R,+) — (L(E),o) 


t = etu 


est un homomorphisme de groupes dérivable sur R, et f'(t) = ue!" = etuu, 


Preuve : 
fG+h)-f() et+hu _ etu cu e*t_ Id ei" (hu)" 
h hk hk hk si n! 
n—1,n h°— 2 uw? 
PANE L _ _- eu + eh NT 
n>1 1e n>2 


L'expression e'"h tend vers 0 quand À tend vers 0, et la fonction qui a À 
n— u? 


associe 32,0 À —— est bornée quand |h| < 1 (on peut en effet écrire la 


majoration ||} ,>2 Du I < > LP dès que |h| < 1). Par conséquent 


lim <E+h) FE) tu 
h—0 L 
Théorème 3 Soient M : R — L(E) etx:R— E deux fonctions dérivables 
(resp. de classe C1). Alors l'application 
g: R — E 
t + Mit)x(t) 


est dérivable (resp. de classe C?) et g'(t) = M'(t)x(t) + M (t)x' (t). 


— € U. 


Preuve : La démonstration est calquée sur celle de la dérivabilité d’un produit 
(qu'il faut savoir retrouver). Pour t fixé, il s’agit de montrer que 


A = st tn 0 — (M'(t)z(t) + M (t)x' (t)) 
tend vers 0 quand À tend vers 0. On a 
À = LUN + h) (+) ME] -(M' (Or (+ M(b x (t) 
_ M (+ n EXD E0 | NON NO, (5 


= M(t+h) as s'(e) 


M(t+h)- M (t) 
+ [Se 


— M! 0] z(t)+(M(t+h)-M(t))x (+). 
Majorons |[A|| par la somme des trois normes, et donnons-nous un réel £ € R*. 
La fonction k + |[|[M (t+ h)]|| est bornée au voisinage de 0 car M est continue 
en t. Il existe donc a € R* et 7, € R* tels que |M(t+h)|| < a pour tout À 
vérifiant ||h|| < 7,. Par hypothèses, on a 


æ(t+h)-x(t 
DR EETTES 
dès que |] < 
LE +R) — M (9 
h 


dès que |A] < n3, et |M(t+h) —- M(t)| < € dès que |h] < 74. Finalement, 
si l’on pose 7 = Min (71,72, 13,4), on obtient bien l'implication 


UOIES 


Ve>0 >0 [Al<n=1IAI<(a+lx(61}+1z (0) es 


2 Systèmes linéaires d’ordre 1 sans second membre 


2.1 Théorème fondamental 


Un système différentiel linéaire d’ordre 1 sans second membre est la 
donnée d’une équation (H) : x! (t) = ux (t) où u € L(E) est fixé. L’inconnue 
x (t) est une fonction dérivable { + x (t) de R dans Æ. Toute fonction déri- 
vable {+ x (t) de R dans Æ vérifiant (A) est appelée solution du système 


différentiel (4). Si l’on rapporte E à une base e = (e1,…,e) et si l’on note 
(&;j) la matrice de u dans e, le système différentiel (H) s'écrit 


z'(t) = Mzx(t) 


où æ(t) = t(xi(t),.…,x9 (t)) désigne le vecteur-colonne des coordonnées de 
x (t) dans la base e. In extenso, le système (A) s'écrit 


din (€) = D ja An5Tj (6). 


Le problème de Cauchy relatif à l'équation différentielle (H) et aux con- 
ditions initiales (t0,x0o) € À x ÆE est celui de la recherche des solutions du 
système (1) vérifiant x (to) = «o. 


Théorème 4 

1) Le système (H) : x'(t) = ux(t) admet la fonction x (t) = elt-to)u (x5) pour 
unique solution satisfaisant la condition initiale x (to) = xo. 

2) L'ensemble S des solutions de (H) est un K-espace vectoriel de dimen- 
sion n. 


Preuve : 1) La fonction x (#4) = elt-to)u (x5) est clairement solution de (H) 
et vérifie x (to) = to. Si x (t) est une autre solution de (A), il suffit de poser 
g(t) = e ‘tx(t) et de calculer g'(t) = —ue ‘tx(t) +e-t“x/(t) = 0 pour 
pouvoir affirmer que g (t) est une fonction constante. En posant g(t)=ceE, 
on obtient x (t) = e!“ (c). La condition initiale x (to) = xo nous donne la valeur 
de la constante : c = e-{04 (x0). Finalement 


æ(t) = er (Cas (xo)) — et—to)u (to) . 
2) L'ensemble S est clairement un K-espace vectoriel et l'application 


wo: EE — S$S 


zo + æ(t)=e"( 


T0) 


est un isomorphisme d’espaces vectoriels. En effet, & est linéaire, surjective 
puisque toute solution x dans S' s'écrit x (t) = el!" (x0) avec ro = x (0) d’après 
le 1), et injective car 


VER zx(t)=e"(xo) =0 = x9 =0. 


2.2 Structure des solutions complexes 


Soit u € L(E). Notons A, …, À, les valeurs propres de uw, puis @1, …, 
a» les multipicités correspondantes de ces valeurs propres dans le polynôme 
ca-ractéristique y, de u, et enfin N(;) — Ker(u— \;1d)" le sous-espace 
caractéristique associé à À. On sait que (voir [1]) 

(1) dim N (À) — di; 

(2) E=N(M)8.….8N(h), 

(3) N() est stable par u. 


Tout vecteur x de E se décompose en x = y1 +. + y avec y; € N (À), donc 


p P 
ett (x) _ > € (yi) _ > et\ildét(u—xld) (yi) 

i=1 i=1 
P œil k 

_ Se #* (u 10 (y) (car etM1d = ex ]d) 
i=1 k=0 
p 

_ pare (t) etñi 
i=1 


œil CTANE 
ot = Y, GX) (y) t* est un polynôme en t à coefficients dans E. 


On peut énoncer 


Théorème 5 (Structure des solutions complexes) 

Les solutions de (H) : x’ (t) = ux(t) sont de la forme x (t) = 3 Ÿ_, P,(t) ef 
où les À; désignent les valeurs propres de u et où les P,(t) sont des polynômes 
de EÎt] de degrés deg P; < @;, œ; représentant la multipicité de À; dans le 
polynôme caractéristique de u. 


Théorème 6 Si toutes les valeurs propres de u sont simples, les solutions 
de (H) sont exactement les fonctions x (t) = 5%, yie!\ où y; décrit le sous- 
espace propre Ker (u — À;1d) de u associé à À. 


2.3 Calcul pratique des solutions complexes 


On remplace x (t) par l'expression 3, P, (t)el\i dans l'équation (H), pour 
obtenir 


p 
SUP) + XP, (6) ef = Vu (P; (6) el* 
i=1 i=1 
EN(X\) EN(X) 


Les appartenance à N (};) se vérifient en remarquant que tous les coefficients 
des P; sont dans N(;), N (}) étant stable par u et donc par u — À;1d. La 
somme directe E = N (A)@...@ N (A) permet d’écrire 


P;(t)+ XP; (6) =u(P,(#)) 


pour tout indice #, et cela déterminera parfaitement les polynômes P;. Cette 
méthode de résolution est dite par identification”. 


2.4 Solutions réelles 


Si (A) est un système différentiel réel, la matrice de u est à coefficient réels. 
Le polynôme caractéristique de « l’est aussi et ses racines sont soit réelles, soit 
conjuguées. Notons A1, …, À, les valeurs propres réelles de u, @1, …, a, leurs 
multiplicités, 4, …, 44, Hi, …, 7, les valeurs propres complexes non réelles et 
Vis 7 les multipicités de w1, …, u,. L'application u : R° — R? est aussi 
une application de C” dans C”, et l’on connaît parfaitement les solutions 
complexes de (H). Pour obtenir toutes les solution réelles de (H), le plus 
simple est d'utiliser le 


Lemme 1 La fonction x (t) est une solution réelle de (H) : x’ (t) = ux(t) si 
et seulement si c’est la partie réelle d’une solution complexe. 


Preuve : Toute solution réelle est bien la partie réelle d’elle-même. Ré- 
ciproquement, toute solution complexe de (A) s'écrit x (t) = y(t) +iz(t) et 
y! (t) +iz! (t) = uy(t) + iuz (t) entraîne y (t) = uy(t), autrement dit y(#) est 
encore solution de (H). 

Une autre méthode possible consiste à donner à priori la forme des solutions, 
puis à ré-injecter l'expression dans l'équation différentielle. Cette méthode, 
dite par identification”, est basée sur le Théorème suivant : 


Théorème 7 (Structure des solutions réelles) 
Les solutions réelles de l’équation (H) : x’ (t) = ux (t) sont de la forme 


= >. P, (t)et\i + ÿ Qi (t) ets + > QOer 
il i=1 


i=1 


où P,(t) € R'[t] sont des polynômes de degrés deg P; < ai, et Q;(t) € Cf] 
avec deg Q; < 7;. 


Preuve : Tout vecteur x € E se décompose de façon unique en 


Da D de un Le (1) 
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avec y € N(X), x € Nu) et z € N(%;) donnés par le Théorème 5. Si 
z € R?, 
Creatas bem 7 (2) 
Les équations (1), (2) et le Lemme ci-dessous entraînent y; = 7; et 2; = 7! 
pour tout 2, soit y, € R et 
T=Y +... + Yp + A +. + 2 +21 +. + Za. 
Toute solution réelle s’écrira donc 
p q q 


et (x) = Net (us) + Def (2) + Set (3). 
1 


i=1 i=1 


En reconduisant le même calcul qu’au Théorème 5, on trouve 


P q 


q 
z (t) — > A0 et + ÿ Qi (t) elhi + Y, R; (t) ei 
j i=1 i=1 


avec 


de la forme désirée. 1 


Lemme 2 On a N{u) = N(f%) pour toute valeur propre 11 de u. 
Preuve : Les racines y et 7 ont même multiplicité a dans v,,, donc 


zeN(u) & (u—uld)" (x) =0 & (u— Hd) ()=0 & TeN(p), 


ce qui entraîne N (4) C N (x). Réciproquement, tout y € N (f) s'écrit y =T 


et les équivalence ci-dessus entraînent x € N (x), soit ye N(u). 


2.5 Etude des systèmes différentiels 2 x 2 


Etudions les trajectoires des solutions de l’équation 


(H) x (t) = Mx(t) 


où la matrice M = ( . ; ) est à coefficients réels. 


a) Si M possède deux valeurs propres distinctes À et u (À < y), alors M est 
diagonalisable et les solutions de (1) sont les fonctions x (t) = ce “uj-+coeltus 
où (uj,u2) désigne une base de vecteurs propres de M. La trajectoire des 
solutions est l’arc paramétré 


{ a (t) = ae 


x (t) = coelt 


et l’on obtient l’une des allures de courbes suivantes : 


FL 


€ = 0, & > 0 et u <0 a =0,@>0etu>0 CC £0et À <u<0 


Si c1c2 £ 0, l’origine est stable si À < y < 0 (cela signifie que l’on est amené 
vers l’origine lorque le paramètre { croît), et instable si 0 < À < y. Dans les 
deux cas, les trajectoires sont les mêmes mais sont parcourues en sens inverse. 


b) Si M possède deux valeurs propres réelles confondues, on envisage deux 
Cas : 
e Si M est diagonalisable, les solutions de (H) sont les fonctions 


— — 
= ce tu + ce “uà 


SI : : 
où (ui,u2) est une base de vecteurs propres de M. La trajectoire est une 
demi-droite de vecteur directeur cui + cous. 


e Si M n’est pas diagonalisable, alors dim Ker (M — Ad) = 1 et 
E = Ker (M — Ad) 5 Ker(M — Ad) =E (À). 


Si (ui, u2) est une base de E telle que uj € E (À) et u2 € E\E(X), une solution 
de (A) s'écrit 
x(t) = e"(xo) 
— gAtet(u-AMd) (cyu1 + cou) = eau + ec (us + t(u — Id) u2) 


Comme (u — Ad)us € E(\), il existe a € R tel que (u — \d)u2 = aur, et 
l'on obtient x (4) = (c1 + coat) e'uj + e*cu2. Les trajectoires des solutions 
admettent la paramétrisation 


À>0 À<0 


c) Si M possède deux valeurs propres complexes non réelles À = a +18 et 
À = a — i, la forme générale des solutions réelles est x (4) = yet + 7e" où 
y = a +ib, a,b € R? et y € E (A). Comme 7 € E (À) et C?=E(N@E (À), 
le couple (y, 7) est une base de C? et 


(0) = Rey) = (RIT) 


sera une base de C? formée de deux vecteurs réels. Dans cette base (a, b) : 
x (t) = (a+ib)eteit + (a — ib)e%e Pt = 2e% (a cos Bt — bsin Bt) 


et la représentation paramétrique d’une trajectoire t + x (t) sera de la forme 


Le diagramme est stable si à < 0, instable si & > 0, et la trajectoire est une 
ellipse si à = 0. 


3 Systèmes linéaires d’ordre 1 avec second membre 


Donnons-nous une fonction f : R — FE continue. Un système différentiel 
du premier degré avec second membre est une équation différentielle de 
la forme 


(E) z'(t)=ur(t) + f(). 
L’équation sans second membre associée est 


(H) z'(t)=ux(t). 


Elle est encore appelée équation homogène associée. 


3.1 Résultats fondamentaux 


Théorème 8 (Méthode de variation des constantes de Lagrange) 
Il existe au moins une solution de l'équation (E) de la forme x (t) = e“t (c(t)) 
où c(t) est une fonction à déterminer. 


Preuve : On a x’(t) = ue" (c(t)) +e“t (c!(t)) et 
(E) & (CH) =fE & CO =e "(f(). 


La fonction { + e-“{(f(t)) est continue sur tout R, et il suffit de prendre 
c(t) = de (f (t)) dt pour conclure. 


Théorème 9 Toutes les solutions de (E) sont obtenues en additionnant une 


solution quelconque de l'équation sans second membre (H) à une solution par- 
ticulière de l'équation avec second membre. 
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Preuve : Si y(t) est une solution particulière de (E), 
{ v'(t) = uxr(t) + f (4) 
y CE) = uy(t) + f () 


— 
& (c-y) (#9 =u(x—y)() 
& x—ysolution de (E) 


z'(t)=ux(t) + f(E) 


& +eE x=y+e"*(c).m (x) 


Corollaire 1 (Unicité du problème de Cauchy) 
Il existe une et une seule solution de (E) satisfaisant la condition initiale 
æ (to) = to. C’est x (t) = y(t) + e"U-to) (x9 — y (to)). 


Preuve : Les solutions de (Æ) sont les fonctions æ(t) = y + e“!(c) où 
c € E et où y(t) désigne une solution particulière de (E). On détermine 
le vecteur c pour avoir æ(to) = xo. On trouve xo = y(to) + e“# (c) d’où 
c=e "0 (xo—y(to)). m 


Le Théorème 9 montre à quel point la recherche d’une solution particulière de 
l’équation (Æ) est importante. Dans la pratique, la recherche d’une solution 
particulière de (Æ) se fait soit par la méthode de variation des constantes de 
Lagrange (Théorème 8), soit par la méthode d'identification exposée dans la 
Section suivante. 


Remarquons bien que la recherche d’une solution particulière de 


a (4) = ux (4) + f(E) + g(t) 


peut s’effectuer en superposant plusieurs solutions particulières d’autres équa- 
tions plus simples. Cela signifie qu’on cherche d’abord une solution particulière 
y(t) de x’ (t) = ux (t) + f (t), puis une solution particulière z (t) de l'équation 
x’ (t) = ux(t) + g(t), pour finalement constater que y(t) + z(t) est solution 
particulière de x’ (t) = ux (t) + f (t)+g (€). Un exemple complet de résolution 
est donné à la Section 4. 


3.2 Second membre de la forme e“Q (t) 
Théorème 10 Soient pu EC et Q(t) Ee Et]. L’équation différentielle 
(E) 2'(t) = ux (t) +e#Q(t) 


admet au moins une solution particulière de la forme e“tP (t) où P est un 
polynôme de degré deg P — deg Q+mult (4), l’entier mult (u) désignant l’ordre 
de multiplicité de u dans le polynôme caractéristique de u. 
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Preuve : Notons r = mult(u), Q(t) = 5% qtf et P(t) = DE prtk où 
px et gx désignent des vecteurs de E. La fonction x (t) = el* P (t) est solution 
de (E) si, et seulement si, 


pet? P (t) + et P!(t) = ue“tP (+) + elQ (+) 
ce qui équivaut successivement à 


(u— uld) P(t) = P'(t) -Q() 


m+r m 
Du nd) (pe)t* = DIE +1) peu — gxlt” 
k=0 k=0 


Vk € {0, l,...,m +r} (u = u1d) (px) — (k “ 1) pr41 — k-: (x) 


e Sir —0, l’endomorphisme u — 114 est inversible et (+) équivaut à 


Pm = (u — HId) ” (am) 
Pm—1 — (u — uld)” (MPm — nt) 


Dans ce cas, le polynôme P est unique. 


e Sir > 0, décomposons p4 et q, dans la somme directe E = N(\)®...@ 
N (À) des sous-espaces caractéristiques : 


{ Pk = PkA + + DR 
k = Ah +. + x: 


Supposons 4 = À, et notons uw; la restriction de u à N'(X;). L’équation (+) 
équivaut à 


Vi VkRE{0,1,.,m+r} (ui —uld) (px) = (k+1)prri,x — x. (+) 


Si à Z 1, l’'endomorphisme u; —u1d est inversible et le système (++) possède 
une solution d’après le premier cas. Si à = 1, alors (u — u1d)" = 0 et il suffit 
de choisir po,, quelconque et définir les p4,1, = px, tels que 


1 
Vk € {0, ls m} Pk+ti,u — k+1 [Qui — uId) (Pk,u) Su ul 
pour que les m + 1 premières équations de (+x*) soient satisfaites. Ensuite, 


m+1.1 étant imposé, il suffit de voir que l’on doit seulement vérifier 
Pm+1,u pose, q 


1 


Vkefm+l,.,m+r} Prin = a 


(ui — Hd) (pk) 
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C’est possible. On remarque surtout que 


1 
Pm+r+1,u — ( 


Gr Done nes (4 Hd) Gmtiu) = 0 


puisque (ui — u1d)" = 0, ce qui prouve que le nombre de coefficients py,u est 
fini et achève la preuve. On a ici trouvé une infinité de polynômes P répondant 
à la question. 5 


Remarque : Le Théorème 10 peut être utilisé dans les cas où le second 
membre f (t) est de la forme Q (t) cost, Q (t)sint, Q (t) cht ou encore Q (t) sht. 


Il suffit de superposer les solutions en notant que cost — te etc. 


4 Exemple 


Cherchons les solutions complexes du système différentiel avec second membre 


x! = 5x — 3y — 3z 
(E) 4 Y =-3r+5y+3z+e 
2 = 9x — Jy — 7z. 


On doit d’abord chercher les solutions générales du système sans second mem- 
bre 
x = 5x — 3y — 3z 
(A) 4 y = —-3x + 5y + 3z 
2 = 9x — Jy — 7z. 


Le système (H) s'écrit X’ = MX où 
5 —3 -3 “ x (t) 
M — —3 5 3 et X=| y = | y{ 
9 —9 —7 z 2 (£ 
Le polynôme caractéristique de M est 


= X 5 3 
Mel = bex à St 
Qu = 


Toute les solutions de (1) seront de la forme X (t) = Ae-t + (B + Ct)e?t où 
À, B, C' sont des vecteurs-colonnes de C°. Pour déterminer les vecteurs À, B, 
C', on injecte cette forme générale des solutions dans l'équation différentielle 
pour obtenir 


AGP ECer +0 (BLCDe" 
—Ae"t + (984 0)e" 4 20te 


MAe * + (MB + MCt)e* 
MAe t + MBe’t + MCte’ 
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La famille des fonctions{t + te} étant libre, on déduit 


meN,aeC? 
MA = -A (M+I)A=0 (M+I)A=0O (1) 
MB=-2B+C &4 (M-21)B=C e{ (Be (2) 
MC =2C (M -21)C =0 (M-21)Ÿ B=0. (3) 


Il s’agit donc de résoudre (1) et (3), puis de déduire C par (2). On trouve 


6. 88 x 0 
G) & | -3 6 3 y |={0 
G 0 6 


T —=ûi 
& A=| y | = a où a € C 
z —3a] 
De même 
à =ÿ BS\" fr 0 
(a) + —3 3 3 y = | 0 
9 —9 —9 2 0 
—9 9 9 x 0 
& 9 —9 —9 y = | 0 
—27 27 27 2 0 
z b1 + ba 
Nr — y —= bi où b1,b2 € C 
z bo 
puis 
3 —3 —3 bi + bo 0 
C=— —3 3 3 b: = | 0 
9 —9 —9 bo 0 


En conclusion, les solutions du système (A) seront 


dé + (b1 + b2) e?t 
X(ti=Ade Fa (Bac aie? + bie?t 
—3a1e ? + be?! 


où a1,b1,b2 € C. Comme on s’y attendait, l’espace des fonctions solutions de 
(H) est de dimension 3 sur C, et une base de cette espace est 


gt e2t e?t 
tr ét ss tm | et |; tn 0 
—3e"t 0 e?t 
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Les solutions générales de (E) s’écrivent comme somme d’une solution parti- 
culière de (Æ) et d’une solution générale de (A). Cherchons donc une solution 
particulière de 


X'=MX+Jet oùJ=| 1 |. (+) 
0 


Comme —1 est racine d’ordre 1 de y (À), on sait que cette solution parti- 
culière peut être recherchée sous la forme X = (A + Bt)e-t. En remplaçant 
dans (+), on trouve 
(B-A-Btje* = M(A+Btje +Je* 
B-A-Bt = MA+J+MBt 


soit 


(M+DA=B-J 
{ (M+1I)B=0. 


Ce système équivaut à 


B=(M+I)A+J (4) 
(M+IŸA=-(M+I)J (5) 
6 -3 -3\°/x 6 -3 -3 0 
(5) & | -3 6 3 y |=| -3 6 3 —1 
9 —9 —6 z 9 —9 —6 0 
18 —9 —9 LA 3 —a — À 
& | —-9 18 9 y |=| -6 | & A- a 
27 —27 —18 # 9 —3a—1 
Prenons par exemple À = {(—1/3,0,—1). Alors 
6 —3 -3 —i 0 1 
B=(M+I)A+J=!| -3 6 3 0 +1 = -1 
9 —9 —6 —1 0 3 


Une solution particulière de (Æ) sera donc donnée par 


(3+8e* 
X=(A+Btje t— —te t 
(—1+3t) et 
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et les solutions générales de (Æ) seront 


ae t+(bi+b)e*+(t- 1e t 
X = ae *t + bie2t — te t 
—3a1e * + bre + (3t— 1)e Ÿ 


On notera que la même méthode peut être utilisée pour résoudre une équation 
dont le second membre est plus compliqué, comme par exemple 


z — 5x — 3y — 3z + cos 2t 
(E) à y = -3r+5y+3z2+et 
# = Ox — Jy — 7z. 


24 0—2t . Le 
Comme cos 2 = EE. le système (F7) s'écrit 


2 


+ 


e-2t 


(a 


X'=MX+Je t+K+L où K = 


(=) ox 


2 
et L — 0 
0 


et une solution particulière de (E') sera obtenue en additionnant les solutions 
particulières de chacune des équations X’ = MX + Jet, X'=MX+Ket 
X'= MX +L. 


5 Equations linéaires d’ordre n à coef. constants 
Il s’agit d'équations de la forme 


CE) 200 (4) + an 1e (4) +. + aox (#) = f (t) 


où les a; sont des constantes complexes et où f (t) est une fonctions continue 
de R dans C. Toute fonction n fois dérivable x (t) de R dans € vérifiant (E) 
est appelée solution de l’équation différentielle (Æ). Le second membre 
de (E) est f (t). L’équation différentielle homogène associée à (E) (en- 
core appelée équation sans second membre) est, par définition, l'équation 


(H) 200 (4) + an_12 D (6) +. + aox (t) = 0. 


Définition 1 L'équation À + an_1X" 1 +. + ao = 0 est appelée équation 
caractéristique de (H). 
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5.1 Utilisation du système différentiel d’ordre 1 associé 


Théorème 11 1) L'ensemble Sx des solutions de (H) est un espace vectoriel 
de dimension n sur R (ou €). 
2) La fonction x (t) est solution de (H) si et seulement si 


où les À; désignent les racines de multiplicité r; de l'équation caractéristique 
X° + an +. + ao = 0, et où ms (t) E R[t] (ou C[t]) est un polynôme de 
degré deg pi < ri. 

3) Le système (tie) D 

4) La solution générale de (E) est obtenue en additionnant une solution par- 
ticulière de (E) et la solution générale de l'équation sans second membre. 

5) Il existe une et une seule solutions de (E) satisfaisant les conditions ini- 
tiales (x (to) , x’ (to) , ne Le best 


est une base de Sy. 


Preuve : 1) x(t) est solution de (AH) si et seulement si 
X (4) = À (oi (6), (6), œn (0) = (x (6,2 (8), 0 (4)) 


est solution du système différentiel d’ordre 1 : 


a (t) 0 1 ” 0 æ1 (t) 
æ9 (4) gi za (t) 
(a - 
0 1 
xt (é) A0 . —An-2 —An—-1 Tn (t) 
D 
X/(0) À X(0) 


L'application 


p: Sx — SH 
at) + KO =t(e(0,2 0,282) 


est un isomorphisme d'espaces vectoriels de de l’espace S7 des solutions de (H) 
sur l’espace S4, des solutions de (H”). Par suite dim Sy = dim Sr = n. 
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2) Le polynôme caractéristique x 4 de À est 


—X 1 0 
—) 1 
À) — 
XA (à) ” 
— À 1 
—4a0 An—2 Gi = À 


= (D (NV + an 1 +. + ao). 


Cela se montre par récurrence en développant le déterminant suivant la pre- 
mière colonne (on applique l’hypothèse récurrente à la seconde ligne) : 


—) 1 0 
—) 1 
xAÙ) = (à)! : * — a (1) 
— 1 
—@1 L nd ni — À 


= (XD A an 2 +. + ai) + ao (—1)" 
(—1)" (x + GATE Hs ao) : 
On sait que la forme générale des solutions de (H) est X (t) = 5?_, P,(t) et 
où P;(t) € ETt] et où deg P; < r;. Compte tenu de l'identification du 1), 
æ(t) aura bien la forme x(t) = 37, p(t)e\t où p;(t) € R[t] (ou C[t]) et 
deg ps < ri. 
3) Le système de n vecteurs (éseit) 


CES = 
1<i<p, 0<s, «engendre l’espace vecto 


riel Sx des solutions qui est de dimension n d’après 1). C’est donc une base 
de S H: 


4) Même démonstration qu’au Théorème 9. 

5) Se donner la condition (x (to) , x’ (to) , ee Coll =, 1) 
revient à imposer la condition initiale X (to) = (1,%2,...,%n-1) à la solution 
de l'équation du premier ordre X’(t) = AX (t) + F'(t), où F(t) désigne le 
vecteur-colonne ? (0,0, ..., f (t)). Le résultat provient de l’unicité du problème 
de Cauchy pour l’équation X’ (4) = AX (t) + F(t) (Théorème 1). 

5.2 Autre méthode 


Si l’on sait que l’ensemble S7 des solutions de l’équation différentielle 


(A) 200 (6) + an_1e D (4) + + aox (#) = 0 
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est un espace vectoriel de dimension n, on peut montrer directement le résultat 
suivant (en conservant les notations de la section précédente) : 


Théorème 12 Une base de l’espace Sn est B — (ésieAit) I née 


Preuve : dimSy = n = #B de sorte qu’il suffise de montrer que B est un 
système libre et que chaque te”! est solution de (H). 

e Montrons que te” est bien une solution de (H). Pour cela, supposons 
que À soit racine de multiplicité r de l'équation caractéristique, et que s < r. 
Soit L le polynôme L (À) = À" + an_1X71 +... + ao. On pose an = 1 et l’on 
définit l'opérateur différentiel 


d _ 
OX 
i=0 
L'équation (H) s'écrit alors L (2) x (t) = 0 et tout revient à montrer le Lemme 


suivant : 


Lemme 3 Sipi(t) est un polynôme et si À € K, alors 


L (£) (p(#) e) = € > 170,0 (4). 


Preuve du Lemme : 


i 


L(R)CO) = Dan (pe) ete 3h04 
i=0 


i=0  k=0 


_: LM po (,) à! = 
nor G-ni 


k=0 i=k 
._ | p (4) d à 
_ Fe . 2 aie (6) ln 
nm k) LrA 14 
D (é) d i 
= AY DC) À 
k=0 k 


n k) 
pi (À 
k=0 


e Le fait que B est un système libre provient du 
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Lemme 4 Si À1, .…, À, sont n complexes distincts deux à deux et si p1 (t), 
.…, Pn (t) sont des polynômes de K [t], alors 


moe... Ln (Dette S pi =.…2= 1, (0) = 0 


Preuve du Lemme : Raisonnons par récurrence sur n. C’est trivial si n = 1. 
Au rang n, l'équation pr (t) et +... + p, (4) et = 0 devient 


pa (#) ef +. + pn1 () eln-1f + pr (E) = 0 


où les coefficients y; = À; — À, sont non nuls et distincts deux à deux. Si 
Pn (t) est nul, il suffit d’appliquer l’hypothèse récurrente pour conclure. Sinon 
posons d = deg p, (t) et dérivons l’égalité d + 1 fois. On obtient 


n—1 


D t)et) 40 2 0 
i=1 
soit 
n—1 d+1 
k Er. 
D Céur (4) pr Fetit = 0. 
i=1 k=0 
L'hypothèse récurrente entraîne ur OC? pe (t) 7 Jim — 0 pour tout 1. 
Si pi (t) n’était pas le polynôme nul, le système (p: (t) , po (os .…, plie Ps) (t)) 


serait libre puisque formé de polynômes de degrés étagés, et l’égalité précédente 
entraînerait pŸT1 = 0, soit s; — 0, ce qui est absurde. On à donc p;(t) = 0 
pour tout à € N,-_1, ce qui entraîne aussi py (t) = 0. La récurrence aboutit. m 
5.3 Second membre de la forme eq (t) 

Théorème 13 Soient q(t) un polynôme de Cft] et y € C. L'équation dif- 
férentielle 


CE) 200 (8) + an12 D (4) + + aox (t) = q (6) el 


admet une unique solution de la forme x + p(t) et où p(t) est un polynôme 
de valuation val (p) = mult (4) et de degré deg p = deg g+mult (4), où mult (y) 
désigne la multiplicité de \ dans le polynôme caractéristique de (E). 


Preuve : Posons d = deg q, s = mult (y) et q(t) = qo + … + qatT. D'après le 
Lemme 3, il s’agit de trouver un polynôme 


d+s 
p(t}= pit +... +parstits = NV prtf 
—=$5 


20 


tel que 


Comme LÜ) (x) = 0 pour tout à < s, cela équivaut à 


n d+s (i) 
LE - pk (k—1)..(k—i+1) ki — q(t) 


n d+s—i ; 
L® Ro | | 
7. Ds. (G+j)(i+j—-1).. (+16 =q(t) 
i=s j—=0 " 
d [inf(d+sin) ,/; 
L®  . | | 
+ >. D... (Gi+j)(G+7—-1)..(j+1)| € =a(t). 


3=0 i=s 


En fait, dans cette dernière somme, ? varie de s à d +5. En effet, p étant de 
degré < d+s, on a p;+; — 0 dès que i +3 > d+s, donc à fortiori si i > d+s. 
Il faut maintenant résoudre 


d+s ; 
LE) 
Vj € {0,...,d} >. 


i=s 


De dede D CobDee. 


il 


Ces d équations forment un système triangulaire dont les inconnues sont les 
Ps, …, pars et dont les coefficients des termes de la diagonale principale sont 
non nuls (en effet L(®) (x) Z 0 puisque y est de multiplicité s dans le polynôme 
L(t)). Ce système est inversible est détermine les ps, …, pq+, de façon unique. 
E 
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